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1 Introduction 


The study of quantum effects in black hole spacetimes comes back to the 
early seventies, when Hawking discovered P] that black holes evaporate by 
emission of thermal radiation (see also W) This result generated enormous 
interest in the subject, especially after Hawking himself speculated |1] that 
the evaporation process will lead to the disappearance of the black hole and 
the information about its formation will be lost forever. This is a radical 
conclusion, as it implies that in the quantum theory the whole process of 
black hole formation and evaporation is nonunitary. 

It is clear, however, that the approximation considered to derive this re¬ 
sult, i.e. the quantization of matter helds in the hxed classical background 
describing the formation of a Schwarzschild black hole, and even the frame¬ 
work used, the semiclassical theory of gravity (see, for instance, u El m), 
cannot lead to a reliable resolution of this paradox. At a certain point dur¬ 
ing the evolution the quantum effects will backreact and modify significantly 
the background geometry, which therefore cannot be considered as fixed, nor 
evolved in a quasi-static approximation. Moreover, once the black hole has 
reached the Planck size quantum gravitational effects will become important 
and cannot be neglected anymore.il Thus it is no wonder that still today 
Hawking’s provocation continues to raise much debate and although most 
of the people do not want to give up unitarity jS] (and between them, re¬ 
markably, now Hawking himself P]), a dehnitive answer on whether and how 
information is recovered in black hole evaporation is still lacking. 

To take into account self-consistently the backreaction effects one needs 
to solve exactly the semiclassical Einstein equations 

( 1 ) 

for the metric Qap, where the quantity on the right hand side represents the 
expectation value of the stress-energy tensor operator of the matter helds in 
a suitable quantum state |T). 

In the hxed Schwarzschild background 

ds^ = = -(1 - rs/rydf + - -- + , ( 2 ) 

_ _ _ (1 - rs/r) 

II To find a way out to this problem one usually argues that quantum gravitational 
effects should always be negligible compared to those due to a large number N of matter 
fields. 
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where = 2GM/c^ is the Schwarzschild radius, three inequivalent quan¬ 
tum vacuum states can be dehned. The hrst is the Boulware state \B) [Til] , 
probably the most natural one, constructed by requiring that in the asymp¬ 
totic region, where the metric becomes minkowskian, it reduces to Minkowski 
ground state |M). It has the property that {B\T^i,\B) vanishes asymptoti¬ 
cally, but the drawback is a strong divergence at the horizon r = rs HH. One 
can circumvent this difficulty by introducing a new quantum vacuum state 
\H), called the Hartle-Hawking state [12], such that is regular at 

the horizon. However one pays a price, i.e. this stress tensor is nonvanishing 
at large r and describes thermal radiation at the Hawking temperature 


T 


hc^ 1 
Gks SttM 


( 3 ) 


where /cs is Boltzman’s constant. The associated physical situation is that 
of a black hole in a cavity, in thermal equilibrium with its own radiation. 
The third possibility, the Unruh state \U) [TSj, is constructed in such a way 
as to reproduce the late time behaviour of the quantum matter helds in the 
classical background of a collapsing star forming a black hole. By requiring 
that no particles are present at inhnity in the past and that {U\Tf^i^\U) is 
regular at the future horizon one hnds, as a consequence, an outgoing flux of 
thermal radiation (the Hawking flux) in the asymptotic future. 

By inserting the hxed background expression |T) in the right 

hand side of Eqs. (CD one can solve perturbatively the backreaction equations 
at 0{h) to hnd the hrst order quantum corrections Sga/s to the Schwarzschild 
metric, i.e. ga/s = gag^ + This is a good approximation to the full 

solution of CD only when 6gag -C i.e. when the quantum terms are 

small compared to the background.** In the case of large mass being the 
Hawking temperature very small this condition is satished, for an evaporat¬ 
ing black hole, for most of its lifetime, but eventually at the late stages of the 
evaporation one faces the problems mentioned above. Turning to the static 
conhgurations, this is a good approximation in the Hartle-Hawking case since 
{H\T^i,\H) never gets large, while for Boulware it is certainly valid at inhnity 
but not at the horizon due to the divergence of the quantum stress tensor 
there. The problem of understanding how, in this case, the quantum ehects 

**The natural coordinates to be used in this context are the Schwarzschild ones given 
in ©■ In more general terms, this approximation scheme is valid whenever the quantum 
terms generate a correction to the curvature which is small compared to its background 
value. 
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modify the structure of the classical horizon is not an easy one. One usually 
disregards this question since Boulware state describes the vacuum polariza¬ 
tion around a static star whose radius is bigger than rs, and therefore this 
divergence is not physically relevant. 

The motivation for our work comes from braneworld physics, where much 
work is being done on the search for black hole solutions in the Randall- 
Sundrum model RS2 HH. This is technically a very involved situation and 
nobody so far has achieved the goal of hnding a hve-dimensional solution 
describing a black hole localized on the brane ng. An interesting physical 
interpretation of this fact comes from the application of the holographic 
AdS/CFT correspondence jT^], for which classical 5D bulk solutions are dual 
to 4D self-consistent semiclassical conhgurations where gravity is coupled to 
quantum matter helds ini. If this is so then a classical static 5D “ black 
hole” will be mapped to a 4D static solution of the backreaction equations. 
Staticity naturally selects the Hartle-Hawking and Boulware states. If, in 
addition, we require that the 4D conhguration be asymptotically flat then 
the choice must be the Boulware state. That this is the correct choice is 
supported by the fact that for large radius the corrections to the Newtonian 
potential in this state match those calculated classically in 5D jl8j . 

In order to solve the backreaction equations in the Boulware vacuum 
the exact expression of {B\T^^{gai 3 )\B) for an arbitrary geometry is needed. 
No such expression exists in four dimensions (an analytic approximation 
for static spherically symmetric spacetimes has however been developed in 
PI). The situation greatly improves if one restricts to the s-wave approx¬ 
imation for spherically symmetric backgrounds. Within this context, in 
section 2 we review the classical and semiclassical theory of gravity cou¬ 
pled to a massless and minimally coupled scalar held. Since in this case 
the expression of {B\Tf^^{gai 3 )\B) is available, we end the section by writing 
down the relevant backreaction equations. These equations cannot be ex¬ 
actly solved analytically, and we approach the problem in two steps. First, 
in section 3, we consider the equations arising in the Polyakov theory, which 
can be derived from ours using a near-horizon approximation for the scalar 
held (this amounts to neglecting backscattering ehects in the propagation 
of the matter helds). These equations are decoupled, in the sense that 
one can derive an equation relating the conformal factor of the metric p 
{ds^ = + dx^) -|- r'^dVL^) as a function of the radius r only (or 0 

or z through the dehnitions r = roe~^ = tqz according to the convenience) 
which is then integrated numerically and from which one can derive the de- 
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pendence on the spatial coordinate p = p{x) and r = r{x). The quantum 
corrected solution is very similar to the Schwarzschild solution from inhnity 
till very close to the classical horizon r = rs, where, as expected, big differ¬ 
ences emerge. In particular, there exists a timelike surface r = tb where the 
radial function r bounces (i.e. the two-spheres reach a minimum radius tb 
and then they start to increase) and, beyond it, a null curvature singularity 
with inhnite radius at a hnite affine distance. Armed with these techniques 
and results we face, in section 4, the full backreaction equations in the s-wave 
approximation, which are much more complicated. The differences with re¬ 
spect to the previous case are that the bounce is located closer to the classical 
horizon r = rs and that the curvature singularity is now timelike and has 
hnite radius. Finally, in section 5 we summarize our conclusions. 

2 Gravity coupled to a massless scalar field 
in the s-wave approximation 

It is convenient, in the context of the s-wave approximation, to work with 
spherically reduced theories. Under the spherically symmetric ansatz 

(is^4) = ds‘f^2) + (4) 

the Hilbert-Einstein action 



reduces to 

= ['•'■R + 2(1 + IVrP)] . (6) 

where the geometrical quantities refer to the radial part ds^ 2 ) of four¬ 
dimensional metric. Note that the radial variable r here plays the role of 
a scalar held with a non-trivial coupling to the radial sector of the metric. 
Einstein’s equations for spherically symmetric conhgurations in vacuum can 
be rewritten as 

2 A S' 

7^# ^ [-2rV,V,r + 9.q2rnr - 1 + |V,f)] = 0 , (7) 

( 8 ) 
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The solution of these equations is the Schwarzschild geometry 


(iS(2) 
r* 

where r* is the so-called “tortoise” coordinate. 

Turning to the matter sector, let us consider the action for a minimally 
coupled massless scalar held (in Gaussian units) 

■S™ = -^/. (10) 


= -(1 


ZUM 

c^r 


) {dt^ — dr 


*2 \ 


2GM, , 2GM, 


( 9 ) 


In the background = gabdx°'dx^ -|- r‘^dVL^ the held / can be expanded in 
spherical harmonics, of which we pick up only the s-wave component 

/ = /(x“) ^ ^Foo- (11) 

r 

Under this assumption, integration of the angular variables in ca leads to 

Sm = ~ j d^x^/^r‘^{V ff . (12) 


Varying this action with respect to the radial part of the metric we obtain a 
two-dimensional stress-energy tensor 


2c 5Sm 


= Z 


ab 


(13) 


which is related to the radial components of the corresponding four-dimensional 
one by the relation 


ab 


Tab 


(14) 


Moreover, by varying oa with respect to r we get the expression for the 
angular components of the four-dimensional stress-energy tensor 


-^ee 


sin^d 


rc 5Sm 
Stta/— 


(15) 
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2.1 Semiclassical theory 

The advantage of the approximation considered is that in this case, unlike 
the full four-dimensional treatment, one can provide an analytic expression 
for the expectation values of all components of the stress-energy tensor. We 
shall briefly review the main steps involved. The details can be found in [7] . 
To this end it is very convenient to parameterize the radial part of the four 
dimensional metric in conformal gauge as 

ds^2) = —e^^dx'^dx~, (16) 

and moreover it is also useful to parameterize the radial coordinate as follows 

r = roe~^. (17) 

One can univocally provide an expression for (T^±) tiy impos¬ 

ing two simple conditions: 

• the covariant conservation laws 

v'‘(rW) = 0, (18) 

which can be rewritten as 

(19) 

• at an arbitrary point X of the spacetime manifold the expectation 
values of the quantum stress-energy tensor (T±±(x^(X))) reduce to 
the normal ordering ones (: T±±{x^{X)) :) when using a locally inertial 
frame based on that point 

{T±±{Cx{X))) = (: T±±(ei(X)) :). (20) 

These two conditions are strong enough to provide a generic expression 
for the expectation values of the stress-energy tensor. In particular, the 
breaking of the classical Weyl symmetry (meaning that classically g°'^Tab = 
0), produces a non-vanishing trace anomaly which can be derived from the 
above conditions. One easily obtains that 

(T) = X(R - e{V4,r + 60,1,) . (21) 
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The full expression for the expectation values of the stress-energy tensor 
components, in an arbitrary conformal coordinate system, is 


{>P|r±±(a:±)|>l.) 


-^id±pd±p - dip) + ^{d±pd±4> + p{d±4>)'^) 
(^|:T±±(a:±):|vl/), (22) 

-^{d+d-p + 3d+(l)d_(l) - 3a+a_0) , (23) 

IZTT 

(^l^l^)p=o - ^{d+d-p + d+pd-(j) + d-pd+(j) 
2pd+d_(j)) . (24) 


The dependence on the quantum state is all contained in the three func¬ 
tions (T| ; T±± : jT) and (\l/|^^|T)p=o- These functions are not independent 
and verify the following relations 

a^(vl/| :T±± : |x[/) + a±0(T|^|vI/)p=o-^9±(a+0a_0-a+a_(/.) = 0. (25) 


2.2 Backreaction equations in the Boulware state 

In a dynamical scenario such as black hole evaporation it is highly nontrivial 
to unravel the precise form of the state-dependent functions (4/| : T±± : 
|T) and (\I'|^^|\h)p=o- However, in this paper we are interested in static 
conhgurations, for which p and 0 are functions of the spatial coordinate 
X = (a;+ — x~)/2 only, i.e. p = p{x) and 0 = 0(a;). This coordinate x 
reduces, in the classical limit, to the tortoise coordinate r* given in ©• For 
the Boulware state it is natural to impose that 


{B|:r±±((,U:|-B) = 0. 


(26) 


This allows to determine, from Eqs. (j25|l . the function 




(27) 


where the index x means derivative with respect to the coordinate x. Thus we 
have all the ingredients we need to write down the backreaction equations in 
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the Boulware state, which describe how the Schwarzschild solution is modified 
due to pure vacuum polarization effects 


where 






SSo 




5(j) 


2c 6S„ 


V^^9 
2c 6S, 


±± 


= (^|T±±|vI/) , 


V^S9' 




±± 


6Sq SSm I T \ 


„3„2 -2(/) 

° - 2d±pd±(j) - {d±(j)Y) 


G 

„3„2 -2(/- 
^ '0^ 

G 


4rQ 


„3„2„-2</. 

— 'o^ 

G 


{d+d-p + d+(j)d-(j) — d+d-(j)) . 


(28) 


(29) 


It is convenient to fix the constant scale tq as follows: = a/A = ■ 


^ Planck _ 

1277 ~ 


m 


127rc^ * 

The static differential equations corresponding to Eqs. (EHl) can then be 
written as 


4^XX ‘^px4^x 6 [pxx Px ff ^Px^x T 6 p(|/> 2 .j 

g2((/i+p) 

XX — + 


A 


= [pxx - ^(j)xx - 


b — (b^ — o = 

''xx Yx pxx ^ 


3p 

XX H” ^Px4^x ~t“ ^p4^XX ~t“ 


3 (0XX 0a;)a 


(30) 

(31) 

.(32) 


To solve these equations we have to add boundary conditions which, in the 
present context, are naturally given by imposing that for very large r the 
solution approaches the classical one ©, he.. 


P = 2 ^^^^ 


2GM. 


c^r 


= r + 


2GM 


ln(l — 


2GM^ 

c^r ‘ 


r = X 


(33) 

(34) 























We shall investigate how the relations p = p{r) and r = r{x) are modihed by 
(static) quantum effects. 

It is convenient to analyze hrst what happens in a simplihed context, 
dehned by neglecting the coupling of the scalar held with the radial function 
r in the classical matter action m- 

3 Polyakov theory’s approximation 

In this section we shall study a simplihed version of the problem outlined in 
the previous section. As already mentioned we shall replace the action (jl2|l 
by a new one obtained by hxing the radial function r = rs = constant. We 
then obtain 

Smatter = j P • (35) 

This approximation is usually motivated by arguing that, in the vicinity of 
the classical horizon r ~ r 5 = 2GM/c‘^, the wave equation for the scalar held 

(-^ + - nrMrf ) = 0 . (36) 

where V (r) is the s-wave potential 

nO = (l--)^, (37) 

reduces to the two-dimensional free wave equation 

+ = (38) 

This latter equation can indeed be derived by varying the action (jnsi- 

The expression for (T|Taft|T) can be derived in a number of diherent ways. 
Following the arguments of subsection 2.1 one arrives at 

{>I>|r±±|>I-) = + (39) 

{4'|T+_|$) = -^9+S_p, (40) 

which are obtained from dll, dH by neglecting the terms depending on 
0 (note that in this approximation = 0). Note that the above 
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expressions can be also obtained from the effective Polyakov action 


Sp = 


h 

967r 


(fx^/—gR^ . 


(41) 


The dependence on the qnantnm state is contained in the fnnctions (T| ; 
T±±{x^) : |\k), which are taken to be zero in the Bonlware vacnnm, i.e. 


{B\ : Tt±(x^) : |B) = 0 . 


(42) 


In the Schwarzschild backgronnd (dSSl) and dSl) x^ denote the Eddington- 
Finkelstein coordinates and the components of the qnantnm stress tensor 
read 


{B\T^±\B) 

{B\T^-\B) 


247r ^ 2r3 8 


h 


247r 


[1-^] 


rs, rs 
2r3 


(43) 


At inhnity {B\Tab\B) —>■ 0 (where \B) rednces to the Minkowski gronnd 
state \M)), while on the horizon {B\T±±\B) —»• —hc^/7687rM^G^. These 
qnantities are strongly divergent when expressed in Krnskal coordinates U ~ 
e”* /2r’s^ Y ^x+i 2 rs Aguiar on the fntnre and past horizons {H~^ and 
respectively) 


{B\Tuu\B) 


{B\T..\B) 

(r-rs^ 


{B\Tvv\B) 


{B\T++\B) 

(r - rs)2 


(44) 


This, in tnrn, means that qnantnm backreaction effects are strong at the 
classical horizon r = rs- 

The semiclassical eqnations in the Bonlware vacnnm can then be written 
as follows 


—^- (dlcj) - 2d±pd±(l) - ( 9 ± 0 ) 2 ) 

Ji„2p—2(j> 1 

— {-d+d.ct) + 2d+ct)d.<i) + ^ 6 ^^'’+'^^) 

„3„2 -2(/) 

2 —— {d+d-P + 5 + 05-0 - 5 + 5 _ 0 ) 

(jr 


- 5±p) (45) 

0 (47) 
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and repeating the steps that led to Eqs. (EDI), (ED) and E2), we rewrite them 
as 


^XX 4^x ‘^Px4^x ^ ^ i^pxx Ptc) 5 

(48) 

^2(4>+p) 


(pxx 2(p^ + ^ — e Pxx , 

(49) 

4^XX 4^x pxx 0 . 

(50) 


We note that when the right hand side of the above eqnations vanishes, while 
keeping hnite the qnotient e^'^/A = we recover the classical eqnations and 
therefore the Schwarzschild solution. Due to the divergent behavior in (I44|l . 
when r approaches the classical horizon rs we expect non-trivial corrections 
to the classical metric. 


3.1 Decoupling the semiclassical equations 

We shall exploit the fact that Eqs. EH),® and EH do not have terms that 
depend explicitly on the variable x and, also, that (HH|l and EH are homo¬ 
geneous differential equations of order two. This allows to write a decoupled 
equation for the function p(0). We use the relations 


Px P0X ) 

pxx T P^xx 1 


(51) 

(52) 


where the dot indicates derivative with respect to 0. Equation can be 
rewritten as 

_ 1+ 2p-l-e2'^[p-p2] 


1 - 

Moreover, subtracting Eq. EH from we get 

. 4>xx p(2 - 

Equations EH and El allow to obtain the desired equation 

.. p[l-2p + e^^{l-p +p'^)] 

p = -- 

or, going back to the radial coordinate r = 

{l + rpr)pr{2r + \pr) 


Prr 


r2 — A 


(53) 


(54) 


(55) 


(56) 
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The equation for 0(a;) can be derived by combining Eqs. (j49|l and (El 


,2 = _^_ 

^ A(l-2p+e2'/'p2) ’ 


and, equivalently, for r{x) 

^ 2 _ 

dx 1 + 2rpr + Xpr 

In the classical limit (A = 0) Eqs. (1^1 and (IKH|l become 


( 57 ) 


(58) 


2(l+r/3,)pr 

Prr 

r 


and 

^2 

dx 


The general solution to Eq. (El is 


g2p 

1 + 2rpr 


Pc 


1 

2 


ln(Al+-) 

r 


(59) 

(60) 


(61) 


where A and B are two integration constants. The Schwarzschild metric 
can be easily recovered by setting A = 1 (i.e., the metric is asymptotically 
Minkowskian, p —>■ 0 as r —+cx)) and B = —2GM/c^. Integration of (jbOj) 
leads to the identihcation of x with the tortoise coordinate r*. 

A natural thing would be to try to solve Eq. (El perturbatively in A. 
This gives a good approximation to the full solution when the quantum terms 
are small compared to the classical ones. At 0(A) this is true for large r where 


Xpr 



XGM/c'^r^ 

1 - 2GMld^r 


(62) 


but not in the near-horizon region where the quantum terms instead dom¬ 
inate. In this region backreaction effects are strong and cannot be treated 
perturbatively. 
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3.2 Numerical solution 


The differential equation dsni) cannot be solved analytically. It must be stud¬ 
ied numerically and for this it is convenient to rescale the radial coordinate 
and introduce the dimensionless parameter ^ = t/a/A. We then get 

. + + (g3, 

Z"' — 1 


This equation allows to analyze in a non-perturbative way the exact function 
p{z). However, for reasons that will be clear in a moment, the function 
p = p{z) is not single-valued. Therefore we have to study, instead, the 
function z = z{p) which verifies the differential equation 


Zpp 


{z + Zp){l + 2zZp) 
z^ — 1 


(64) 


Imposing as boundary condition that the solution behaves, for very large 
z, as the classical one 

p(z ^ oo) = ^\n{l - (65) 

with a = GM/c^y/X, we can generate numerically the solution to (jEH)- We 
find that the quantum corrected solution is everywhere similar to the classical 
one, up to the vicinity of the classical horizon. We can observe this behavior 
in Fig.l. We have chosen a black hole of small-size (a = 10^) since in this 
case the differences between the classical and the semiclassical solutions can 
be better appreciated. 

We observe that for regions far away from the classical horizon [z » 1, 
i.e. p —> 0) the numerical solution and the classical one are very similar. 
However, in the vicinity of the classical horizon the quantum corrected solu¬ 
tion suffers a bounce, absent in the classical solution, around p ~ —8.3 and 
then grows up slowly. This is the reason why we have had to solve numeri¬ 
cally 2 ; = z{p) instead of p = p{z). Note that this point appears at a finite 
value of p, and, therefore, that the time-time component of the metric gu 
does not vanish. The existence of this bouncing point is better represented 
in Fig.2, where we plot the derivative Zp in terms of p. The existence of a 
zero for Zp signals a bouncing point for the radial function. 

This qualitative behavior of the radial function is maintained irrespective 
of the size of the black hole. For a Solar-mass black hole (a = 10^®) the bounce 


13 





p 

Figure 1: Classical (dashed line) and numerical (solid line) plots of the func¬ 
tion z{p) for a = 10^. 
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appears at p{zb) ~ —91. This result implies, as we will see later with more 
detail, that the classical horizon is eliminated by the quantum corrections. 
However, it is important to remark that the value of the conformal factor of 
the metric at the bouncing surface is very small. In fact the redshift for a 
signal emitted by a static observer at the bounce and received at inhnity 

( 66 ) 

^Bounce 

transforms, for instance, a Planckian energy Esounce ~ into an 

energy of the order ~ 10~^^eV at inhnity. Moreover we hnd that, for 
a = 10^®, Zb = zs + 8.48 x , where zs = 2a = 2 x 10^®. This shows 
that the bouncing surface is indeed very close to the classical horizon (the 
difference between zb and zs increases as one reduces the mass). 


3.3 Behavior of the metric around the bounce 

Around the bounce zb the function z{p) behaves 

z{p) ^ ZB + ^A{p- pb)^ + ... . 

Plugging this expansion into we get immediately that 

Zb 


Therefore 


A = 


p(z) ZZ p(zb) ± 


-B 


I2 {zI-1){z-zb) 
Zb 


and 


Pz 


±i 




2zb{z - Zb) 
To estimate the form of the metric at r ~ 


(67) 


( 68 ) 


(69) 


(70) 


( 71 ) 


we need to hnd the relation between r and x. From Eq. (EH) we have that 


= 1 + 2zp^ + pI = 1 + 2rpr + Xpl 


(72) 
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In the region r ^ tb the right hand side is dominated by the pure quantum 
term = Xp^, with p^ given in ((TO)), and so we have 



(73) 


Therefore the form of the metric is approximated by 



(74) 


The quantum corrected geometry is not singular at the bounce, as it can be 
checked that all curvature invariants are regular at r = r^. Note that r is 
not the good spatial coordinate to extend the metric (|7H) beyond one 
should rather use x (see (USD). 

Finally we should briefly note that the surface r = is not an event 
horizon since gui^B) 7 ^ 0. However, there grri’^s) = 0 and this means that 
< 0 for points x inside tb and only for r = tb we have d^r‘^\rg = 0. 
This means that the surface r = tb still plays the role of an apparent horizon 
for outgoing radiation. 

3.4 The geometry beyond the bounce 

We shall now investigate the geometry beyond the bouncing surface. To this 
end we study the 2D curvature for p < p{zb)- Starting from the expression 
of the 2D curvature 


R = 8e ‘^^d+d-p = -2e , 


(75) 


where pxx is given by (see Eq. (IS2D where here and in the next formulas p 
and p are written in terms of derivatives with respect to z) 


Pxx T Z Pzz^4^x ZPz(f)xx 1 


(76) 


and from (jSHD we get 


R = -2z^ 


zp,^ + {‘2z + p^)pI ^_^p^^ 


(77) 
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We can get a simplified expression for the curvature taking into account our 
basic differential equation (El for p{z) 


R 


-2z‘^ 


_Pz^ 2p 
1 + zp^ 


2z^ 


p^{2z + p^) 
z‘^ — 1 


e 





Moreover from El we have 


\z‘^{l + 2zpz+ pI) \z‘^{l + 2zzp +z'j) 
and therefore the hnal expression for the curvature is 

_ 2 pz{2z + pz) _ 2 (2zZp + l) 

A(^2 _ 1) (1 + 2zpz + pz"^) A(2;2 _ 1) (1 + 2zzp + z^) 


(78) 


(79) 


(80) 


The natural singularity at 2 ; = 1 (r = a/A), which mimics the classical 
singularity at r = 0, does not belong to the physical spacetime because of 
the existence of the bounce encountered before. Singularities can only be 
generated by zeroes of the denominator in dHl. Before we have Zp > 
so 1 + 2zzp + Zp can never be zero. At z = zb we have Zp = 0, so 


R 


2 

A(4 -1) ■ 


(81) 


After the bounce Zp < 0 and, therefore, one can potentially encounter a cur¬ 
vature singularity. The numerical analysis indicates that such a singularity 
can be found only when 2 ; —*• -fcxD. For this we need that (1 -h 2zzp + Zp) 0 
when ^ -(-cxo. The zeroes of the above second order polynomial Zp = —z± 

Vz‘^ — 1 are an exact solution to the differential equation (El. Since the nu¬ 
merical analysis shows that, for large 2 ; in the interior region, Zp ~ —l/2z + ... 
this means that asymptotically we have Zp —z + z"^ — 1 + f{z), where 

z"'f{z) ^ 0 as 2 ; —A -1-cxo for every positive integer number n.b The behavior 
of the curvature is then 


R ~ 


1 

TfW ’ 


(82) 


note that the solution Zp = —z + — 1 corresponds to the one found in m 

as an exact solution to Eqs. (ESll . lO and CT where the exponential term is 

neglected. 
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which clearly shows the existence of a singularity at z = +cx). A detailed 
computation, using (ESI), shows that (for ^ > +oo) 

~ \/^ . (83) 

Using Eq. dZl it is easy to realize that in the limit z > +cxo we have 
f(z) ~ and this implies that the scalar curvature goes to —oo there as 

R~-—. (84) 

Moreover, such a singularity is null (i.e. ^ 0 as z ^ +oo) and is 

located at a hnite affine distance from any hnite Finally we remark that 
the singularity arises due to the the branching point for the radial function 
r = zV^, with respect to the spatial coordinate x, displayed in (jHSI)- This 
is the underlying reason for the generation of the curvature singularity at 
X = —oo. 


4 Quantum corrections in the s wave approx¬ 
imation. 


The approximation used in the previous section consists, essentially, in ne¬ 
glecting the effects of the potential barrier for the wave equation (jnni)- In 
this way we have simplihed considerably the technical problem. It is natural 
at this point to ask whether the results obtained are maintained when the 
effects of the potential are included. The detailed analysis presented before 
has allowed to introduce all the techniques that we shall use to attack the 
full problem in the s-wave approximation. The equations to solve are now 
more involved. However, since the conceptual line to follow should be clear 
we will focus only on the most important points. 

In the Schwarzschild spacetime the expectation values of the stress tensor 
components in the Boulware state (see (E2I)-(Eni) and (EZD) are 


{■BlUilS) = A| 


2rc 15 rL 


fi 


l*> ^ 


247r'' 8 167rr^ 

^ M _ rs-, rs 
127r^ r ^2r3 ’ 


(l-f)'ln(l-f), 


h 


7h rs h r| 

1671- 7-3 27r 


(l-^)(l-^)ln(l-^) . (85) 

lyt 7 ^ 
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All these quantities vanish asymptotically, while on the horizon the leading 
divergence is the same as in the Polyakov case (I44|l . Again, this means 
that to inspect the near-horizon region r ~ we need to solve exactly the 
backreaction equations. 


4.1 Decoupling the semiclassical equations 

We shall now proceed in parallel to Section 3 to decouple the system of 
differential equations (|Sni)-(E2D to generate a single equation for p = p{r). 
Using the relations p^ = p4>x and p^x = 'p{(l>xY + P<Pxx we transform them 
into 


0XX - - ‘^P4>1 

(P{4>+p) 

(pxx — 20^ H-— 

(j^XX 03, P0a; P0XX 


[p0x + p(t)xx - P^ffi + 600^ + 6p02] (86) 

[p0x + P<pxx - 3(0xx - (pl)] (87) 

e^‘^[300^ -f Spcpxx + Qp(pl + 6p0xx 
3 {4>xx — 4>x)x -^ 

2 (px 


From (IHF)|l and (IHZl) we get 

0a;a;(l ~ P^^^) = [l + 20 + e^'^(0 — 0^ + 60 + 6p)] 0 

„2((/)+p) 

0xx(l - (0 - 3)e^'^) = [2 + (0 + 3)e^'^] (pl - 


A 


which can be rewritten as 


(1 - pe^<P)e^(^+p) 


XD 


[l + 20 + e2'^(0 - 02 + 60 + 6p)] e2('^+^) 


XD 


where 


(89) 

(90) 


(91) 

(92) 


D = - [1 - (p - 3)e^*] [1 + 2p + - p + 6f} + 6p)] 

+ [1 - pe^'l'] [2 + (p + 3)6^*] . 

Taking into account that 



(93) 

(94) 
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and plugging the expressions (EH) and EH into EHI) we obtain, after a 
straightforward but very long calculation, a third-order differential equation 
relating p and r 


-I- 3A — 6Ap -I- rpr{2r‘^ -|- 6A -|- rXpr) 


X 


72A^p^(r^ -|- 3A -|- rXpr) -|- 


-|-6rAp^l2rA — pr(6r^ -|- 40r^A -|- 63A^ -|- 2rApr(4r^ -|- 14A -|- rXpr)) + 

+rX{2r‘^ -I- 15A -|- 2rXpr)prr^ + r^^rA(6r^ -|- 56r^A -|- 93A^)p^ -|- 

+2r^X\r^ + 6A)p^ + p^Ar^ + 54r^A + lQ2r^X‘^ + 162A^ - Qr^A^rr) + 
-hr^(2(r^ -h 5r^A -h 12A^)prr + 9X^pl^ — SrX'^prrr) + rpr{Ar‘^ — Qr‘^X — 


12X^ + A(2r^ - 14r2A - 69A2)p^^ 



(95) 


The second factor in the above equation is the relevant one since it leads 
directly, when A = 0, to the classical equation EH- The quantum corrections 
to the Schwarzschild metric should be then computed by exactly solving 
the above differential equation. Introducing the dimensionless coordinate 
2 ; = e“‘^ = r/\fX we get an ordinary differential equation for p = p{z). 
However, as we have already explained in the analysis of the Polyakov theory, 
since the function p = p[z) could not be one-to-one it is more appropriate 
to work directly with the differential equation for the function 2 ; = z{p). It 
reads as follows 


18p(-21 Ap)zz^p -1- 216p^2;^ A- 2z^zl{2zp - Zpp) + 2z'^Zp{3zp -f 2Zp 

-Zpp) + 3z^Zp(^(31 - Ap)zp - 8(3 -h 10p)Zp -F (23 - 4p)zppj 
+%z^zl ((27 - 23p)zp + 12p(l + p)zl - IbpZp^ + 2z^Zp{l + 21 z^ 
—bZpZpp) — 3z^{^ — 3zpp + 2zp{—2 + (—27 -|- 8p)Zp + 2(2 -|- p)zpZpp) 
-Zpp^ + ( 562 ;^ - 6(1 6p)2;J - 9zIp + Zp{lAzpp + 3zppp)^ = 0 . (96) 


4.2 Numerical solution 

We solve numerically the above equation by imposing that, for very large 
2 ; and p —>• 0, the solution approaches the classical one. We hnd that the 
solution is almost identical to the classical one up to the vicinity of a surface. 
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which we also denote as zb, located very close to the classical horizon. The 
result (for a Solar-mass black hole a = GM/c^\/X = 10^®) is depicted in 
Fig. H which shows the existence of a bounce for the radial function at 


2 -10 


38 


1 . 5-10 


38 


5-10 


1-10 


38 


39 


- 90.5 



P 


- 5-10 


39 


Figure 3: Plots of the function Zp{p), classical (dashed line) and numerical 
(solid line) for a = 10^®. 

Pb ~ —89.69, where 



(97) 


and A is a positive coefficient to be computed numerically. Moreover, we 
hnd that zb ^ Zs + 1.76 x 10“^^. Therefore, its location is much closer to the 
classical horizon than in the Polyakov theory approximation of section 3. 

4.3 Branching point for the radial function 

The main difference with respect to the Polyakov theory appears in the re¬ 
lation between r and x. In terms of the function z{p) we have 



2 


2 


(1 -F ZZp) 


(1 -F ZZp) 
Xzz^pD 


zZpD 


(98) 


(99) 
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where 


D = -^{ - [z'^Zp + z + 3zp] [z^z^ - 2z^zl - z^Zpp - bzzj - z^Zp + 6pzJ] 

+ \z‘^Zp + z^^lz^z^p + — z^Zpp + zz^pW . ( 100 ) 

In Fig. (ED we show the typical behavior of the function D in terms of p. 


D 



Figure 4: For a = 10^® the bounce of the radial function {D = +cx)) is located 
at p ~ —89.69. At p ~ —89.72 the function D vanishes. 


Note that in the vicinity of the bounce Zp —>■ 0, and from (cnnD we have 
D Zp^. Therefore 

~ 4 . (101) 

dr Zp 

and hence 

dr 

In other words, the expansion of r in terms of x has to be of the form 
r K, rB + a{x — XBY + ■■■■! where a is a numerical constant, in agreement with 
the analytic behavior encountered in dZSl). 

A similar argument allows to determine the behavior of r in terms of x 
around the zero of the function D. This happens at p = pm ~ Ps, just after 
the bounce. Around the zero of D we have D{r) rM — r. Therefore 

-\/rM - r , 


dx 

dr 
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(103) 








and hence 


r ^tm - /3{x - (104) 

where /? is a nnmerical positive constant. The radial fnnction has a branching 
point at a: = xm, which turns out to be the minimum possible value for the 
“tortoise” coordinate x. 

The form of the metric in this region is 

+ dx"^) + {tm + B{x — XMY^^YdVt^ (105) 

where, according to our previous analysis, the function p(r) is hnite and 
regular aX tm- The above metric has a singularity at r = tm, which cannot 
be avoided by a change of coordinates. It is indeed a curvature singularity 
as we now show. The AD scalar curvature can be expressed, in terms of p 
and 0, as follows 


/?(4) = _2e-2p _ 20,, + 302] + 2e2<> . (106) 

The hrst term is just the two-dimensional scalar curvature R = —2e~'^^pxx, 
which according to dzni), (EH) and dnH is 


2 —z^Zpp + 2z^Zp + z"^ + Qzzp — 6pZp 

Xz'^D zz^ 


(107) 


At a: = xm, where D{xm) = 0, R and also the second and third terms in 
(ITO are divergent. Finally we note, from (HDH, that this singularity is 
timelike, has hnite radius tm and is located at a hnite affine distance away. 


5 Conclusions 

The existence of a bounce for the radial function r, prior to the emergence of a 
spacetime singularity, is perhaps the most signihcative result of our analysis. 
It already appears in the simplihed Polyakov theory and it is still there in 
the most accurate s-wave approximation. It is natural to expect it to persist 
in a full treatment of the problem. 

Our results are perhaps not surprising from the semiclassical point of 
view, where due to the strong divergence of the Boulware stress tensor at the 
Schwarzschild horizon important deviations from the classical behavior are 
indeed expected to arise. However they constitute an important prediction 
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for the search of static 5D braneworld configurations with asymptotically flat 
boundary conditions on the brane: such solutions are not black holes, but 
rather naked singularities. In the 4D semiclassical context this reinforces 
the idea that the Boulware state describes the vacuum polarization around 
a static star, not a black hole. Indeed the natural thing for a black hole is to 
be time dependent and to evaporate via the Hawking effect. 

Finally we point out that our analysis does not exclude the existence of 
static braneworld black holes, but the price to pay is to give up asymptotic 
flatness on the brane. This means, in the dual AD theory, to replace Boul¬ 
ware with the Hartle-Hawking state. 
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